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CN . Abstract 

Wc establish existence and uniqueness results for the singular initial value prob- 
er , lem associated with a class of quasilinear, symmetric hyperbolic, partial differential 
equations of Fuchsian type in several space dimensions. This is an extension of 
earlier work by the authors for the same problem in one space dimension. 
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^ 1 Introduction 

^j. In previous work [1], we have established existence and uniqueness results for the sin- 

gular initial value problem associated with quasilinear systems of partial differential 
equations (PDE) of Fuchsian type, with non-analytic coefficients and in one space di- 
mension. In the present paper, we generalize these results to systems in several space 
dimensions. While much of the earlier work with Fuchsian systems involves PDE sys- 
| terns with analytic coefficients, here (and in [1]) we allow coefficients which are C°°, are 

of finite differentiability, or are elements of certain weighted Sobolev spaces. We recall 
that Fuchsian techniques have found many applications in the study of nonlinear wave 
equations [15], in mathematical physics [16], and in mathematical cosmology and the 
study of Einstein's field equations of general relativity [1,3-10,13,16-18,20]. 

The earliest applications of Fuchsian techniques to mathematical cosmology con- 
^ \ cerned equations with analytic coefficients. One such theorem by Kichenassamy and 

Rendall [17] has been applied to several families of cosmological solutions of Einstein's 
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equations. Existence theorems in the class of smooth (but not analytic) functions have 
proven more difficult to establish. Rendall in [18] relies on energy techniques for sym- 
metric hyperbolic systems and establishes that a sequence of analytic solutions converges 
to a smooth solution; he then applies his theorem to deal with the Gowdy vacuum equa- 
tions with T 3 spatial topology. In subsequent work, this technique has been applied by 
Stahl [20] with partial success in attempts to handle the Gowdy equations with S 2 x S 1 
and 5 3 spatial topology, and has also been applied to a family of polarized T 2 -symmetric 
solutions to the Einstein equations by Clausen [12]. 

More recently, a Fuchsian method has been introduced by Beyer and LeFloch [5-8] 
who define solutions in certain smooth and Sobolev-type function spaces for semilinear 
symmetric hyperbolic systems; they apply their theory to establish the existence of a 
large class of Gowdy spacetimes with T 3 spatial topology and with certain asymptotic 
properties in the neighborhood of the singularity. In their approach, a sequence of regu- 
lar Cauchy problems (obtained by working away from the singularity) is introduced and 
it is shown that their solutions converge to a solution of a singular initial value prob- 
lem (defined below), provided that suitable regularity assumptions in weighted Sobolev 
spaces hold. The follow-up theorems in [1] and in the present paper provide an exten- 
sion of [5-8] to the class of quasilinear equations. Other results in the smooth class are 
found in the literature and require various structural hypotheses of the Fuchsian systems 
under consideration and the space of solutions; cf. [18] for further references and, for the 
Sobolev class, [14,16]. A particular distinguishing feature of our approach in [1,5,6], as 
well as in this paper, is that an approximation scheme is at the core of the method. This 
scheme can be implemented for numerical computations straightforwardly and contains 
useful built-in convergence and error estimates [2,5,7]. 

As mentioned above, the present paper generalizes our results in [1] to an arbitrary 
number of space dimensions. In order to achieve this result, we work with dimension- 
dependent Sobolev embeddings and inequalities, and we track all dimension-dependent 
arguments through the construction. Indeed, Sobolev inequalities are crucial in several 
places in the argument, and we refer the reader to the beginning of Section 3, below, for 
further details. 

The paper is organized as follows. In Section 2, we introduce the class of PDE 
systems of primary interest, referred to as quasilinear symmetric hyperbolic Fuchsian 
systems, and we discuss some basic terminology relevant to the singular initial value 
problem for such systems. We then formulate our main existence and uniqueness the- 
orem; cf. Theorem 2.4. We also include, as an example application, a discussion of 
a nonlinear Euler-Poisson-Darboux type equation; this example serves to exhibit im- 
portant features arising with the Einstein system. Next, Section 3 is devoted to the 
proof of Theorem 2.4 and, in particular, discusses some of the energy estimates needed. 
Concluding remarks are found in Section 4. 
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2 Statement of the main result 



2.1 Quasilinear Fuchsian systems 

We are concerned with the following system of PDEs for the unknown u : (0, 5] x T n — > 

S°(t,x,u)Du + S j (t,x,u)tdjU + N(t,x,u)u = f(t,x,u), (2.1) 

where each of the maps S° and is a symmetric d x d matrix- valued function of the 
spacetime coordinates 1 (t,x) and of the unknown u (but not of the derivatives of u), 
while / = f(t,x,u) is a prescribed M rf -valued function, and N = N(t,x,u) is a matrix- 
valued function of (t, x, u). It is convenient to scale the time derivative operator and set 
D := tdt = ijj = x -^, while dj := for j = l,...,n. For notational convenience, 
we shall often leave out the arguments (t,x), and instead use the short-hand notation 
S°(u), Si(u), N(u),f(u). Notationally, S°(u) (for instance) may be regarded as a map 
u i y S°(u) between two function spaces (as discussed below) and, therefore, may be 
written as S°(u)(t,x). The notation 

L(u)[v](t,x) := S a (t,x,u)td a v + N(t,x,u)v 

is also used throughout this paper, so that the left-hand side of Eq. (2.1) can be written 
as L(u)[u]. Here, a = 0,1, ...n, so that the above summation includes the time derivative 
term as well as the space derivative terms. 

Presuming (as stated explicitly below in Definition 2.2) that the matrices S a and 
N are non-singular in a neighborhood of t = 0, we observe that with our choice of 
coordinates, the PDE system (2.1) is singular exactly at t = in the sense that the 
PDE coefficients vanish there (and nowhere else nearby). Hence the singularities which 
we consider in this paper are "spacelike". Apart from this, the convention that the 
equations are singular at t = does not result in a loss of generality, since for solutions 
whose singularity locus is rather at t' = tp(x), one could always introduce the new time 
coordinate t = t' — tp(x). 

At this point the reader may wonder why the zero-order term N(u)u is included in 
the principal part and not in the source f(u). We leave these terms separate since, by 
convention, / includes terms of "higher order" while the term N(u)u contains terms of 
the same order as the principal part in t (as t — > 0). Specific conditions on N(u) and 
f(u) are stated and discussed below. 

In order to measure the regularity and the decay of both the solutions and the 
coefficients near the singularity t = 0, we introduce a family of weighted Sobolev spaces. 
Letting [i : T n — > M. d be a smooth 2 function, we define the matrix 

Tl\j\{t,x] := diagfV^), ...,t~w( x )) (2.2) 

1 While all results here are stated for functions defined on the spatial manifold T n , our techniques 
and results could be extended in similar form (with suitably adjusted function spaces) for other spatial 
manifolds, including R™. 

2 By "smooth", we mean that sufficiently many derivatives exist and are continuous, ft is straightfor- 
ward to check how many derivatives are necessary in each argument and this is left to the reader. 
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and use it to define the following norm for functions w : (0,6] x T n i— > R d : 




SUp ||^[/i]«;||H8(T») 

te(o,<5] 



(2.3) 
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here H q (T n ) denotes the usual Sobolev space of order q on the n-torus T n , a denotes 
a partial derivative multi-index, and the standard Lebesgue measure is used for the in- 
tegration. Note that this norm only controls spatial derivatives. Also note that the 
derivatives d a operate on both the matrix TZ\p] and w; if fi is not constant, then loga- 
rithms in t to the power \a\ have to be controlled. 

Next, we define the function space Xs^, q (T n ) to be the completion of the set of 
functions w 6 C°° ((0, 6] x T n ) for which the above norm is finite. A closed ball of radius 
r about in Xs^ iq (T n ) is denoted by Bs^ q ^(T n ). Note that in our discussion below, 
we include the argument T n only if we wish to emphasize the dimensional dependence; 
we shall often write X$^ q in place of Xs tf i, q (T n ), with the argument understood to be 
T n . To handle functions which are infinitely differentiable and for which we control all 
derivatives, we also define the space X^oo '■= D^=o ^ n ^ ne following, we refer to 

quantities /i as exponent vectors and we always assume that these are smooth. We 
write v > \jl for two exponent vectors (of the same dimension) if each component of v is 
larger than the corresponding component of (i at each spatial point. 

In working with d x d coefficient matrix- valued functions (such as S a ), we use anal- 
ogous norms and function spaces. In these cases, each component of the normed matrix 
is weighted by t raised to the negative of the corresponding component of a smooth 
exponent matrix £ : T n — > M. dxd (analogous to the exponent vector \x above). For 
convenience, we use the notation Xs£ q (T n ) to denote these function spaces. 

2.2 The singular initial value problem 

In contrast to the Cauchy problem for Eq. (2.1), which seeks a function u that satisfies 
Eq. (2.1) and that equals a specified function ur to i : T n — > M. d at t = to > 0, the singular 
initial value problem seeks a solution of Eq. (2.1) with prescribed asymptotic behavior 
in a neighborhood of t = 0. More specifically, one prescribes a "leading order term" uo, 
which may be either a function or a formal power series on (0,6] x T n , and one looks 
to find a solution u such that w := u — uq suitably decays relative to u$ at a prescribed 
rate in a neighborhood of t = 0. Substituting u = uq + w into Eq. (2.1), one obtains a 
PDE system for w, which takes the form 



where 3~(uo)[^] := /(uq + w) — L(uq + w)\uq]. The problem of existence and uniqueness 
for the singular initial value problem is now equivalent to establishing the existence and 
uniqueness of a solution w to Eq. (2.4) with a prescribed asymptotic behavior as t — > 0. 



L(uq + w)[w] 



F(uo)M 



(2.4) 



4 



Using the spaces Xg „„ to prescribe the asymptotic fall-off rate, we now state the singular 
initial value problem in the following form. 

Definition 2.1 (The singular initial value problem). For a given choice of a leading 
order term uq and the parameters 5, u, q, the singular initial value problem consists of 
finding a solution w £ Xg„ q to the system Eq. (2.4) and, therefore, a solution u = uq + w 
to the system Eq. (2.1). 

How does one choose a leading order term? One method, described in [16], is to 
find a power series in t with coefficients depending possibly on spatial variables, such 
that the most singular terms in the equation are canceled. Another method, which has 
proved useful in cosmology and in applications to the Einstein equations, is to select 
a leading order term as an exact solution to a simplified system which is derived from 
the Einstein equations by dropping the spatial derivative terms relative to the time 
derivative terms. This approach leads to the so-called velocity term dominated (VTD) 
system, and using the Fuchsian method we can in certain cases verify the existence of 
solutions with asymptotically velocity term dominated (AVTD) behavior. We refer to [1] 
for a general method for doing this which provides a formal series of leading order terms 
and generalizes easily to n dimensions. We make no further comments on the choice of 
leading order terms in the rest of this paper. 

2.3 Quasilinear symmetric hyperbolic Fuchsian systems 

We establish (cf. Theorem 2.4, below) that, under certain conditions, there exists a 
unique solution to the singular initial value problem as described in Definition 2.1. 

Definition 2.2. Fix positive real numbers 5 and s, non-negative integers qo and q, and 
an exponent vector fi : T n — > Mr, together with a leading-order term uq : (0, 5] x T n — > R d 
(with so far unspecified regularity). The system Eq. (2.1) is a quasilinear symmetric 
hyperbolic Fuchsian system around uq if there exist matrices Sq(uq) that is positive 
definite and symmetric and independent of t, Sq(uq) (j = 1, . . . ,n) that are symmetric 
and independent of t, and Nq(uo) that is independent of t, all belonging to the space 
H qo , and if there exist vector functions fij : T n — > M. d (j = 1, . . . , n) with strictly positive 
components, such that each of the "remainder matrices" 

S?(«o + w) :=S°(u +w)- S° (u ), 

S{(u + w) :=n\Pj - 1] (S> + w)- Sg'(«o)) , 

Ni(u + w) :=N{u + w)- N (u ), 

maps all elements w € Bg> „ q s for all 5' € (0, 5] to elements S®(uo + w),Sl (uo + w) and 
Ni(uq + w) in Bgi Cqn where Q is some exponent matrix with strictly positive entries, 
and r > is some radius. One also requires that Si(uq+w) and S{(uq + w) be symmetric 
matrices for all w £ Bgt n q s . 
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It is important clarify the notation used here and below. While the quantities Sq(uo) 
and iVo(uo) are explicitly time-independent, they do depend on the t — >■ behavior of 
the leading order term uq. For convenience below, if a choice of leading order term has 
been fixed, we may omit the explicit dependence on uq, and write simply Sq and iVo; 
the dependence on the leading order term is then implicit. We use the same notational 
shorthand with Sf (uo + w) and N\(uo + w), omitting the explicit dependence on u$ so 
long as the choice of the leading order term is fixed and unambiguous. We note that 
Sf (uq+w) and Ni(uo+w) are both explicitly time-dependent, with the time dependence 
coming from uq as well as from w. We also note that so long as the system Eq. (2.1) is 
a quasilinear symmetric hyperbolic Fuchsian system as defined above in Definition 2.2, 
the perturbing coefficients 5f,-/Vi belong to the space X^^,q, with Q positive and thus 
are non-singular at t = 0. 

It follows from Definition 2.2 that if Eq. (2.1) is a quasilinear symmetric hyperbolic 
Fuchsian system, then it is also symmetric hyperbolic for all t £ (0, 5] (so long as 5 is 
sufficiently small so that S° is positive definite). Hence for smooth initial data (or for 
data in H q (T n ), with q > n/2 + 1) prescribed at to G (0,<5], the Cauchy problem is 
well-posed in the usual sense (away from t = 0), with solutions belonging to the space 
C(I,H q (T n )) for a sufficiently small interval / C (0,5]; see, for instance, [21]. We note 
that since solutions to the Cauchy problem are only defined for t bounded away from 
the singularity at t = 0, we know nothing a priori regarding the singular behavior of 
these solutions, nor whether they belong to some space X$^ q . 

Our existence and uniqueness results rely on further structural conditions, including 
the following one. 

Definition 2.3 (Block diagonality with respect to fj,). Suppose that uq is a given leading- 
order term and ji is an exponent vector. The system (2.1) is said to be block diagonal 
with respect to ^ if the following commutators vanish (for all w = u — uq E X&^ q ) 

S (u)TZ[fi}-TZ[fi}S (u) = 0, 

s j (u)n[ij} -n[n]s j {u) = o, 

N(u)H\ji] - K[n]N(u) = 0, 

where TZ[fj] is defined in Eq. (2.2), and if the same property holds for all relevant spatial 
derivatives of S° (u) , (u) , N (u) , as well. 

This condition is essential in deriving our energy estimate below in Section 3. It 
ensures that the principal part operator only couples those components of the unknown 
which decay in t at the same rate. Another quantity which plays a role in the derivation 
of energy estimates is the energy dissipation matrix 

M := SqVo) diagGui, .., fi d ) + iV (tt ), (2.5) 

which is a function of the spatial coordinates x, only. 

We can now formulate the central result of this paper, as follows. 
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Theorem 2.4 (Existence and uniqueness for the singular initial value problem). Suppose 
that Eq. (2.1) is a quasilinear symmetric hyperbolic Fuchsian system around uq (with a 
choice of the parameters 5, s, fj,, q and qo as specified in Definition 2.2) and is block 
diagonal with respect to fi. Suppose also that q > n/2 + 2 and qo > n/2 + l + q. Then, for 
some 5 G (0, 5], there exists a unique solution u to Eq. (2.1) whose remainder w := u — uq 
belongs to Xg ^ with Dw G Xj _ v provided the following structural conditions are 
satisfied: 

(i) The energy dissipation matrix Eq. (2.5) is positive definite at all spatial points. 

(ii) The map 3{uq) : w h-> /(uq + w) — L{u$ + u/)[«o] is well-defined, and, for every 
5' G (0,5], it maps B$i^ q ^ s to Xy v q for some exponent vector v > \i. 

(Hi) For all 5' G (0,5], for some constant C > and for all w,w G Bginq, s , one has 

\\3(uo)[w\ - J(uo)[w]\\ 5 '^ q <C\\w- w\\s',n,q, 

and, the constant C below being independent of 5' , 

\\H U o)[ W ] ~ 3(uo)[w]\\6' ,v,q-l + \\Si(u + W) - Si(UQ + w) \\s'X,q-l 
n 

+ W S i( u o + w )~ S i( u o + w)\\ S ',c,q-i + \\Ni(u + w) - iVi(u + ^)||5',c, g -l 
< C\\w - w\\ s >, M -i. 

Furthermore, if the conditions above hold for all integers q > n/2 + 2 and qo > n/2 + l + q, 
then the remainder w belonqs to X? , while Dw G Xr 

Analogous statements can be derived for time derivatives of arbitrary order, by taking 
time derivatives of Eq. (2.1). 

This theorem is mainly a generalization of the corresponding theorem in [1] to ar- 
bitrary spatial dimensions n. However, there is also a small change in Conditions (ii) 
and (iii): While in [1], Condition (ii) must be checked for all w in Xg^q, here this 
condition must hold necessarily only for those w contained in an arbitrarily small closed 
ball Bgngg. A similar change is made here for Condition (iii). This change is important 
in order to encompass operators 3~(uo) for which Condition (ii) is only satisfied if w is 
bounded in some way. As a particular example we consider a singular initial value prob- 
lem, for which ^(uq) : w i— > 1/(1 + w) — (1 — w). If w G Bg jlM q,s for a choice [i > and for 
a sufficiently small s, then one can show that 3"(uo)[u>] G Xg : 2^,q and hence Condition (ii) 
is satisfied. However, Condition (ii) is violated if s is too large, and hence the version of 
the theorem in [1] does not apply. A sketch of the proof of Theorem 2.4 is presented in 
Section 3. 

We observe that in the hypothesis for Theorem 2.4, the degree of regularity required 
for the leading-order matrices Sq(uq), Sq(uq) and No(uo) (specified by qo), and that of 
the perturbation matrices Si(uq + w), S{(uq + w) and N\(uq + w), differ slightly from 
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that of the remainder w (given by q) . The regularity required for the asymptotic data uq 
(implicitly specified in Condition (ii)) is also slightly different. These gaps arise in the 
course of our proof; we note that they arise in particular in working with the higher order 
energy estimates and the corresponding Cauchy problems for derivatives of w which are 
needed to control the regularity of solutions. It is not clear if these gaps can be avoided 
through the use of a different method of proof. In any case, the gaps vanish in the C°° 
class of solutions, corresponding to q = qo = +00. 

In its applications (see, for instance, Section 2.4 below), this theorem often allows 
one to find an open set of values for the exponent vector /i for which there exists a unique 
solution to the singular initial value problem. A lower bound for this set 3 can originate 
in Condition (i), while an upper bound is usually determined by Condition (ii). Both 
bounds on the set of allowed values for \x provide useful information on the problem. 
The upper bound for fi specifies the smallest regularity space and, hence, the most 
precise description of the behavior of w (in the limit t\,0), while the lower bound for 
(j, determines the largest space in which the solution u is guaranteed to be unique. We 
note that this uniqueness property must be interpreted with care: under the conditions 
of our theorem, there is exactly one solution w in the space X^nq-, although we do not 
exclude the possibility that another solution may exist in a larger space, for example, 
in X$^ q with ju < \i. Note that even if a given system does not satisfy our hypothesis 
above, there exists a systematic method which allows one to "improve" a leading-order 
term uo; cf. the discussion of (order-n)-leading-order terms given in [1], or Section 2.2.4 
in [16]. 



2.4 Example 

It is instructive to discuss an example PDE system which has some of the features 
of the systems used to model cosmological solutions to Einstein's vacuum equations. 
We consider the following first-order system in two spatial dimensions (with spatial 
coordinates (x,y)): 

D Ul - u 2 =0, (2.6) 
Du2 — Aii2 — u\td x U3 — uitdyU^ = t e 113 (e Ul — U2U4) , (2-7) 
U1D113 — U1U3 — u\td x U2 =0, (2-8) 
u\Dui — U1U4 — u\td y U2 =0, (2-9) 

where A > and e > are given parameters. This system is derived from the following 
second-order, quasilinear Euler-Poisson-Darboux (EPD)-type equation 

D 2 u - XDu - ut 2 (u xx + u yy ) = tHd x u [e u - t 2 d t ud y u) , (2.10) 

provided we define the new unknowns 

u\ := u, U2 := Du, 113 := td x u, := td y u. (2-11) 



3 A number A is denned to be a lower bound for the allowed values of the vector /x if each component 
v a of v satisfies the condition v a (t,x) > A for all x in the domain of [i. A similar definition holds for an 
upper bound for \x. 
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To recover solutions of Eq. (2.10) from solutions of Eq. (2.6)-(2.9), we must impose 
the third and fourth of the relations in Eq. (2.11) as constraints on the choice of initial 
data (for the Cauchy problem) or constraints on the choice of asymptotic data (for 
the singular initial value problem), and then verify that the evolution preserves these 
constraints. However, we ignore this issue in this discussion. 

Based on considerations discussed in [5], one might guess that a solution of Eq. (2.10) 
is likely to behave as h + t x g near t — > 0, for h and g a pair of M-valued asymptotic 
data functions. The definition of first-order variables Eq. (2.11) (and the notion of 
ODE- leading order terms from [1]) suggests the choice 



u = (h(x,y) + t x g(x,y),Xt x g(x,y),td x h(x,y),td y h(x,y) 



(2.12) 



as a parametrized set of leading order terms. As we determine below, this leading 
order term, for appropriate choices of h and g, together with Eq. (2.6)-(2.9), comprise 
a quasilinear symmetric hyperbolic Fuchsian system (in the sense of Definition 2.2). 

To show this, we substitute u = uq + w into the system Eq. (2.6)-(2.9), expand in 
components of w = (^1,^2,^3,^4), and then comparing with Eq. (2.1) we find that 



S°(u + w) =diag(l,l,/i 



t x g + wi,h + t x g 



wi), 



S x {u + w) 



and 



N(u 



w) 






V 

/ 



V 





Vo 






t x g 







t x g 



\ 

t x g -wi 







w\ 









/ 



-1 

-A 











t x g 



t x g 









t x g 



\ 



J 
\ 

- Wl J 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



where the source term reads 

f(u + w) = (0, t\td x h + w 3 ) (e h+tXg+m - (Xt x g + w 2 ){td y h + w 4 )) , 0, 

The expression of the matrices Sq, 5q, and iVo follows directly by setting t = and 
wi = in the previous expressions. We then determine that f3 x = f3 y = 1. We also see 
that h = h(x,y) must be positive so that 



Sg = diag(l,l,M) 
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is positive definite, and we require that h G H qo (T 2 ) with qo = 3 + q; we require q > 4 
in order to satisfy Theorem 2.4 with n = 2. Then, Sq, SQ,and iVo belong to H qo (T 2 ) (see 
Definition 2.2). 

We seek solutions with remainder w G X$^^ q for an integer g as in Theorem 2.4 and 
for some exponent vector \x = (jUi, /i2, A*3, /m)- in order to interpret u> as the remainder of 
a solution with the leading-order term above, we require that //i,/i2 > A and ^13, ^4 > 1. 
The block diagonal conditions, however, force all components of /i to be the same. We 
therefore label each component simply by the same symbol \x and hence require that 
fi(x,y) > max{l, A} for all spatial points. The energy dissipation matrix Eq. (2.5) is 



M = 5 ° 



ffi -1 \ 

fi- A 

ii- 1 

\0 At-1/ 



'mi 



which is positive definite if /J-(x,y) > max{l, (A + \/l + A 2 )/2} for all spatial points; this 
is a stronger restriction than /i > max{l, A}. 

We now check Condition (ii) of Theorem 2.4, and compute 

y(«o)H 

/ \ 

t e (td x h + W 3 ) (>+* A 9+«-i - (\t X g + iD2)(t^/i + w 4 )) + t 2 (/l + Wi)(/l ra + 

At A+1 (/l + ^l)fe 

v At A+i (^+^i)^ y 

Here we see that if h G H q+2 (T 2 ), g G H q+1 (T 2 ), and 10 G (for g > 4 as above) 

and so long as n(x,y) < min{2, 1 + e, 1 + A} at each spatial point (x,y) G T 2 , we find 
3"(no)[w] G Xs :U ,q for some exponent vector v > fi. 

Let us combine the implications of Condition (i) and Condition (ii) of Theorem 2.4. 
For q > 4, we require that h G H q+3 (T 2 ) be a positive function, and g G H qJrl (T 2 ). The 
upper and lower bounds for the exponent vector fi can be satisfied if and only if 

/15 3 + 8e + 4e 2 \ 
0<A<nnn(-, ^ ). (2.17) 

Finally, Condition (iii) in Theorem 2.4 follows immediately from the following two 
lemmas, which were established in Appendix B of [1] and readily extend to general n > 1. 

Lemma 2.5. Fix any integer q > n/2 and let 4>i,(f>2 : M d — s- M be q-times continuously 
differentiable functions such that w \— > <pi(w) is a map Xg ^ q — > Xg viq and w \— > (f>2( w ) is 
a map Xg^ q — > Xg iU2>q , for (scalar) exponents v\, vi, and exponent vector fx. Moreover, 
suppose that for each s > and for all 5 G (0, 5], there exist constants C±, Ci > (which 
may depend on s but not on 5) such that 

[w] - <t>i[w\\\s, Vl , q < Ci\\w - w\\s, m , 

II02N - Mw]h,u 2 , q < Gt\\w - *ll w 
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for all w,w G Bg „ gs - Aou> ; Zei := <fii4>2- Then it follows that w \— > 9(w) is a map 
Xs,n,q — > ^5,u 1 +u 2 ,q- Moreover, for each s > and for all 5 G (0, 5], there exists a constant 
C > (which in general depend on s, but not on 5) such that, for all w,w G Bg „ ga , 

\\o[w] - e[w\ h iV1+va>q < c\\w - w\\s^ q - 

Lemma 2.6. Suppose q > n/2. Let ipw := expollj, where IT : M, d — > R is the projection 
to the ith component of d-vectors. Then w i-> ip^(w) is a map Xg t ^ q — >• Xg^^. Moreover, 
for each s > and for all 5 € (0, 6], there exists a constant C > (which in general 
depends on s but not on 5) such that, for all w,w G s , 

||^ ) H-# ) NII 5 - 0l9 <c'|k-*|| 5 - Ml? . 

Theorem 2.4 now yields the following result. 

Theorem 2.7 (Existence and uniqueness for our example system). Let h G H q+3 (T 2 ) 
be a positive function and g G H q+1 (T 2 ) with q > 4. The system Eqs. (2.6) -(2.9) 
with parameters e, A > satisfying Eq. (2.17) admits a unique solution u = uq + w 
with uq given by Eq. (2.12) and with w G Xgj^^^^^ for fi(x,y) G (max{l, (A + 
Vl + A 2 )/2},min{2, 1 + e,l + A}) for all (x,y) G T 2 . 

Comparing Theorem 2.4 and Theorem 2.7, we note that while in the former there is 
no direct prescription of the regularity of the leading order term in the hypothesis, in the 
latter the hypothesis requires that h and g (which make up uo, as stated in Eq. (2.12)) 
lie in the specified function spaces. This is consistent, since the regularity requirements 
imposed on S§ and Nq in Theorem 2.4 in fact imply regularity restrictions on uq; we see 
this explicitly in this example. 

A particular problem of Theorem 2.7 is that we only construct solutions with re- 
mainders w in restricted spaces Xg^^ q . In particular, uniqueness only follows in 
a restricted sense: the full space for the remainder which is compatible with the singu- 
lar initial value problem associated with the leading-order term Eq. (2.12) is given by 
Mi) M2 > A, and (13, fj,4 > 1, but here we only obtain uniqueness of solutions in a subspace 
of this. The full uniqueness result can be obtained by means of (order n)-leading order 
terms [1]. 

We remark that if for instance A > 2, it is in fact reasonable to expect that for the 
leading-order term Eq. (2.12) one cannot guarantee the existence of a solution to the 
singular initial value problem. The heuristic reason is as follows. The leading-order 
term Eq. (2.12) is obtained under the assumption that the source-term and the spatial 
derivative terms in the equations can be neglected at t = and that, to leading-order, 
(ill, U2, U3, 7/4) is therefore a solution of the resulting ODE system. However, since for 
example the second spatial derivative terms in Eq. (2.10) are 0(t 2 ), which is of the same 
order as the leading-order term for A > 2, this assumption breaks down. 
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3 Proof of the main theorem 



3.1 Preliminary 

The proof of Theorem 2.4 given in this section closely follows that of the main theorem 
for systems in one spatial dimension which is presented in [1]; therefore, we refer the 
reader to that paper for the details and only focus here on the main steps and point out 
key differences which occur for space dimension n > 1. 

There are two main steps to the existence proof. First, we prove an existence theorem 
for a linear singular initial value problem of Fuchsian form and, second, we use this 
result to formulate a contraction mapping argument for the full quasilinear problem. 
To prove existence to the linear problem we use a sequence of Cauchy problems with 
the initial times tj characterizing each of them approaching t = 0, and we use energy 
estimates to show that the sequence of solutions to these Cauchy problems does converge 
to a solution to the singular problem. We start by defining the relevant class of linear 
Fuchsian systems. To simplify the notation for the linear theory, we set uq = (without 
loss of generality), so that the solution u coincides with the remainder w. 

Which aspects of the proof of Theorem 2.4 are most affected by the replacement 
of one space dimension by more than one? Since Sobolev estimates and embeddings 
depend on dimension, wherever these play a role, dimension matters. This is true for 
the well-posedness theorems for the Cauchy problem for symmetric hyperbolic systems 
(see, for example, [21]), for handling non-smooth coefficients in the linear singular initial 
value problem (see Proposition 3.5), and also for the contraction mapping argument. 
The first of these affects our formulation and derivation of our energy estimates; see 
Lemma 3.3 and Lemma 3.4, below. The second requires us to work with products of 
elements of Sobolev spaces; the resulting dimensional dependence can be seen in the 
Moser estimates (stated for instance in Proposition 3.7 in Chapter 13 of [21]). The use 
of Sobolev estimates in carrying out the contraction mapping argument results in the 
requirement that the regularity index q satisfy the condition q > n/2 + 2 (consistent 
with the requirement that q > 3 for one space dimension). 

3.2 Linear symmetric hyperbolic Fuchsian systems 

As noted, the starting point for our proof of Theorem 2.4 is a study of linear symmetric 
hyperbolic Fuchsian systems. We define these now, obtain energy estimates for them, and 
then use those estimates to prove existence and uniqueness of solutions to the singular 
initial value problem for linear Fuchsian systems. 

Definition 3.1 (Linear symmetric hyperbolic Fuchsian systems). Suppose that 5 and r 
are positive real numbers, q and q$ are non-negative integers, [l : T n — > R d is an exponent 
vector, and ( : T n — > M. dxd is an exponent matrix such that £ > 0. The system Eq. (2.1) 
is a linear symmetric hyperbolic Fuchsian system if the following conditions are 
satisfied: 
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(i) The matrices S° , , and N are independent of u and can be decomposed as 

S°(t,x) =S° (x) + S° 1 (t,x), 

S*(t,x) =K[1 - Pj(x)] (Sfa) + S? (*,*)) , 

N(t,x) =N (x) + N 1 (t,x), 

for matrices Sq(x) (positive definite and symmetric) , Sq(x) (symmetric) and Nq(x) 
in H qo , and for matrices Si (t,x) (symmetric), S^(t,x) (symmetric) and N\(t,x) 
in Bg£ qr . Here /3j : T n — > M. d (j = 1, . . . ,n) is a collection of exponent vectors 
with strictly positive components. 

(ii) The constant 5 is sufficiently small so that S°(t,x) is uniformly positive definite. 
By uniformly positive definite we mean that S°(t,x), given by the decomposi- 
tion in Condition (i), is positive definite uniformly 4 at all t G (0, 5] in the I? -sense 
with respect to x, and for all S® S Bg^ q>r . 

(Hi) The source term is linear in the sense that 

F(u )[w] :=f(t,x,u) = f (t,x)+F 1 (t,x)w (3.1) 

with fo G ^6,v,q an d the. matrix F\ satisfies TZ[^\FiTZ[^\~ l £ Bg^ q ^ r . Here v is a 
exponent vector with v > a. 

Both in the linear and in the quasilinear theory, the matrices Sf and Ni are consid- 
ered to be perturbations of Sq and No. In formulating Theorem 2.4, we seek estimates 
that are independent of the particular choice of Sf and N±. This can be done provided 
these quantities are bounded; hence the introduction of the balls Bg^ q r - Throughout 
the proof we keep careful track of which quantities the constants in various estimates are 
allowed to depend on. To make this precise, it is useful to have the following definition. 

Definition 3.2. Suppose that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system 
for a chosen set of the parameters 5, u, £, q, qo and r. Suppose that a particular estimate 
(e.g., the energy estimate Eq. (3.3)), involving a collection C of constants, holds for 
solutions of Eq. (2.1) under a certain collection of hypotheses H. The constants C are 
defined to be uniform with respect to the system and the estimate so long as the following 
conditions hold: 

(i) For any choice of S®, S{, N±, and F\ in the perturbation space Bg^ >q>r as in Defi- 
nition 3.1 which are compatible with the hypothesis H, the estimate holds for the 
same set of constants C. 

(ii) If the estimate holds for a choice of the constants C for one particular choice of 
5, then for every smaller (positive) choice of 5, the estimate remains true for the 
same choice of C. 

4 Since we assume in addition above that 5o is positive definite and the perturbation Si is bounded, 
this implies that 5° is positive definite pointwise at all (t, x) 6 (0, 8] x T n for all S? £ B$x,q,r if qo and 
q are sufficiently large (thanks to the Sobolev inequalities). 
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As mentioned above, the first step in our proof of Theorem 2.4 involves setting up a 
sequence of Cauchy problems (with the initial times approaching zero), and controlling 
the solutions of these Cauchy problems using energy estimates. In carrying this out, it is 
convenient to make the temporary assumption that the perturbation quantities and 
Ni and the source functions fo and F\ belong to smooth subsets of their respective spaces 
(cf. Definition 3.1). Presuming this, we say that Eq. (2.1) has smooth coefficients. 
A later step in the proof removes this smoothness assumption. We note that although 
this assumption - that Eq. (2.1) has smooth coefficients - implies that the quantities 
Sf,Ni,fo and F\ are differentiable to all orders, it does not not guarantee that all 
derivatives have controlled asymptotic behavior. This control holds only for a set of 
derivatives, as labeled by the relevant function space. 

We now define the energy functionals of interest, state the energy estimates, and 
discuss the proof of these estimates. Generalizing the corresponding notion in [1], we 
define 

E^ n [v]{t) := \e- Ktl (5°(t,.)W(*, 0, WCVMV))^) ( 3 - 2 ) 

for any function v : [t ,5] x T n ->• R d (with v{t, •) G L 2 (T n ) for each t G [t ,5]). In 
the above definition S° is the matrix appearing in (2.1). It is easily proved that these 
energy functionals are equivalent to the L 2 (T n )-norms with equivalence constants which 
are uniform in the sense of Definition 3.2, if S° can be decomposed as Condition (i) of 
Definition 3.1. Based on this definition, we have the following energy estimate. 

Lemma 3.3 (Fundamental energy estimate). Suppose that Eq. (2.1) is a linear sym- 
metric hyperbolic Fuchsian system for the parameters 5, /i, £, q, qo, r (according to Defi- 
nition 3.1), has smooth coefficients, and is block diagonal with respect to [i, with q = 
and qo > n/2 + 1. Suppose also that the energy dissipation matrix Eq. (2.5) is positive 
definite for all x £ T n and, in addition, DSi,diS{ G Bg£Q a for all i,j = 1, ... ,n for 
some s > and some exponent matrix £ with strictly positive entries. Then for any 
initial data vr tQ i G H qo (T n ) specified at some to G (0,(5], there exists a unique solution v 
to the corresponding Cauchy problem, and there exist positive constants k, 7 and C such 
that v satisfies the energy estimate 

y/E^[v](t) < E^ Kri [v](t)\ t=t0 +C f s^WKlriis, -)fo(s, OHl^t*)^ (3.3) 

J to 

for all t G [to, 5]. The constants C, n, and 7 may be chosen to be uniform 5 . In particular, 
if one replaces Ur to i specified at to by any vi tl ] specified at any time t\ G (0, to]; then the 
energy estimate holds for the same constants C , k, 7. 

The proof of Lemma 3.3 is essentially the same as that for the one-dimensional 
case [1]. We note that the existence of unique solutions to the n + 1 dimensional Cauchy 

5 While the constants C, k and 7 here can be chosen to be uniform in the sense of Definition 3.2, 
there generally does not exist a choice which holds for all <5, Sq, Sq, No , f3j ,r, s, f , £, fj, and v. 
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problem corresponding to Eq. (2.1) (which follows from, e.g., Proposition 1.7 in Chap- 
ter 16 of [21]) plays a key role in this result, and that the inequality for qo stated in the 
hypothesis of Lemma 3.3 is needed to guarantee such existence. 

We also need to control higher order spatial derivatives of the solutions, for which 
we establish the following energy estimate. 

Lemma 3.4. Consider a linear symmetric hyperbolic Fuchsian system for the given 
parameters 5, /x, (, q, qo and r, which satisfies all of the conditions in Lemma 3.3, except 
that here we allow for arbitrary integers q > 1 and qo > n/2 + 1 + q. Assume as well 
that 6 DSf € B$^fl :S . Then there exist positive uniform 7 constants C,p such that for all 
e > 0, the solution v(t,x) of Lemma 3.3 satisfies 



- e](VMV)||if9(T«) < C(||ft[/i](t ,-M l|jf«(T") 



+ 



to 



ifa-i(T«), 



ds 



(3.4) 



The same choice of constants C and p can be used for any initial time to £ (0,(5). 

The inequality for go comes again from the condition for well-posedness for the 
Cauchy problem in n spatial dimensions, but now also from the fact that in deriving the 
above estimate we take q spatial derivatives of the coefficients. 

Proof. We only discuss the case q = 1, noting that the argument generalizes to higher 
values of q. Let di be any fixed spatial derivative operator i = 1, ...,n. As in the case 
n = 1, we derive an equation for diV by acting on both sides of Eq. (2.1) with the 
derivative operator di. Doing this for all values of the index i, we obtain the following 
system 

S° • Ddv + tSi ■ tdjdv + N ■ dv = % + F x ■ dv, (3.5) 
for the "unknown" (n • d)-vector dv := (div, . . . , d n v) T , the (n ■ d) x {n ■ cQ-matrices 



5° := 



(S° 
S° 

: 



Vo ... 

the (n ■ d)-vector 

with components 
Jo 



0\ 



5°/ 



Si 









Si 



Si J 



N 



(N 






N 



0\ 



Nj 



70-) 

JO ! 



f{n) 



6 We note that the condition diS{ € Bs,£,o 
choice q > 1 



difo + (PiF! - 8iN)v - (OiS^rHfo + Fi« - iVf), 

n) of Lemma 3.3 is now implied by the 



i,3 = 1 

7 We note however that C and p generally depend on q 
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and the (n • d) x (n ■ d)-matrix 
with components 



( if' j) 



F^ j) ■= FiStj - tdiS j + tdiS ^ )-^^ 

All these terms can be handled in essentially the same way as they are in the n = 1 case 
in [1], thereby producing the inequality Eq. (3.4). □ 

To prove that solutions to the singular initial value problem for these linear systems 
exist, it is useful to work with weak solutions: We first show that weak solutions exist, 
and then prove that the weak solutions are differentiable enough to be strong solutions. 
To define weak solutions for the linear symmetric hyperbolic Fuchsian system Eq. (2.1), 
for 5, /i, C) Q-, Qo an d r, we first define a test function (f> : (0, 5] x T n — > W 1 to be any 
smooth function with the property that (p(t, •) = for all t > T E (0, 5). Then, we define 
the operators C and T acting on w E X$ tlM fi (similar to [1]) by setting 

(C[wU) := - ft ((K[tJ,}S w,Dct>) + ^{KMS^tdrf) - (K[(,}Nw , <f>) 

JO JT n v ~T 



i=i 







and 



+ {n[n\ (s° + ds° + ni^DTi^ny^s 

n 

+ {K[f4~HtdjK[f4)S j + tdjS'Yjw^dxdt, 



(F[w},0):=[ [ (K[ fJ ](f + F 1 w)w,<l))dxdt. 
Jo Jr n 



In [1] we prove a lemma which shows that these operators are well-defined bounded maps 
on Xs jfMj Q-, this lemma easily extends to the n-dimensional case. We call w E ^8,^,0 a 
weak solution of the equation Eq. (2.1) if it satisfies {C[w] — F[w], (ft) = for any test 
function (p. 

Suppose now that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system for 5, 
fi, C, Q, Qo an d r as in Definition 3.1 and suppose that it is block diagonal with respect 
to n and that the energy dissipation matrix is positive definite. Working with q = 
and qo > n/2 + 1, we recall the construction of the approximate solutions from [1,5]: 
Letting {t m } be a monotonically decreasing sequence of times in (0, 5} with limit zero, 
we define v m : (0, 5] x T n — > R d to be the piecewise function which is identically zero 
on (0,t m ] and which satisfies the Cauchy problem for a linear symmetric hyperbolic 
Fuchsian system Definition 3.1 with zero initial data at t m on [t m ,<5]. It is shown in 
detail in [1] for the one spatial dimension case that, under appropriate conditions, this 
sequence of approximate solutions converges to a weak solution of the singular initial 
value problem (as noted above, we presume here that the leading order term for this 
problem vanishes). The same argument works here for general n. The energy estimate 
Lemma 3.3 is crucial to this argument. This construction of weak solutions gives rise to 
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a solution operator EI which maps fo to a weak solution. The arguments above depend 
on the system having smooth coefficients, so in particular, the operator EI to this stage is 
only defined for fo E -^5,^,0 n C°°((0, 5] xT n ). However, since EI is a linear and bounded 
map Xs t u,o H C°°((0, 5] x T n ) — > Xs^ t o, it can be extended uniquely as a linear bounded 
operator to the full space fo E X$ Ui0 . We note that at this stage of the discussion, we 
retain the smoothness assumptions of the other coefficients above. 

The next step of the argument, which is again discussed in detail in [1] for the n = 1 
case, is to show that for q > 1 and qo > n/2 + 1 + q, the energy estimates of Lemma 3.4 
allow one to verify higher regularity of the solution. In particular, the solution operator 
EI can be shown to be a well-defined map X$ u ^ q — > X§ One can also show that the 
weak solutions w are differentiable with respect to t and that Dw E Xs )fjL) q~i- It then 
follows that w is a strong solution of the (linear) system. 

We next proceed to remove the smoothness assumption which has been imposed 
on 5f , N\ and F\. This step is the first one where the generalization to more than 
one spatial dimension requires some care, as we discuss in the proof of the following 
proposition. 

Proposition 3.5 (Existence of solutions to the linear singular initial value problem). 
Suppose that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system for 5, p, C, Q, 
qo and r as in Definition 3.1 and suppose that it is block diagonal with respect to p for 
qo > q + n/2 + 1 and q > n/2 + 1. Suppose that the energy dissipation matrix Eq. (2.5) 
is positive definite. Then for all fo E Xg t uq with v > p there exists a unique solution 
w : (0, 5] x T n — > M. d to the singular initial value problem with zero leading order term 
such that w E X$^^ q and Dw E Xs^ t q-i ■ The solution operator EI : fo i— >■ w satisfies 



for some positive uniform constants C , p. 

We note that this result also holds for the case qo,q = oo: If the conditions of this 
proposition are satisfied for all integers qo > q + n/2+l and q > n/2 + 1, then w E ^5,^,00 
and Dw E X$^ tOQ . However, the g-parametrized sequence of constants C and p occurring 
in the estimate of the solution operator may in general be unbounded as q — > 00. 

Proof. We restrict attention to the case F\ = 0. The basic idea is to set up a se- 
quence of equations [w] = fo defined by a sequence of smooth coefficients S^ 
and ^V"i[i] which converge to the non-smooth coefficients in the Xs^ y q norm. For each 
of the corresponding singular initial value problems, we can apply the (smooth co- 
efficient) arguments discussed above, and we thereby obtain a sequence of solutions 
E Xs t ^.q, Dw\i] E Xs^^^i to this sequence of singular initial value problems, with 

< C" 1 1 /o|| • 

We now argue (as in [1]) that the sequence converges in Xs^^-i- To do this, we 
note that the function — wyj is a solution of the system 




(3.6) 




(3.7) 
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where we consider the right-hand side to be a given source term. Now, using the proper- 
ties of the solution operator corresponding to L^, we determine that the || • H^g-i-norm 
of the left-hand side of Eq. (3.7) can be estimated by the || • H^q-i-norm of the right- 
hand side, for some v > u. To estimate the right-hand side, we use the block diagonal 
condition and the decompositions of the coefficient matrices which follow from Defini- 
tion 3.1 of a linear symmetric hyperbolic Fuchsian system. To estimate the resulting 
products which consequently appear on the right-hand side, we apply the Moser (Sobolev 
product) estimates, which require that q — 1 > n/2. Then, the fact that the sequence 
of coefficient matrices is a Cauchy sequence implies that wu\ is a Cauchy sequence in 
Xs,tM,q-i- We show that the limit w £ X$^ q _i is a weak solution, and in fact a strong 
solution of the equation with Dw S Xs tfJ , t q-2- 

The remaining step for the existence argument is to show that we can recover the 
loss of regularity and in fact verify that w G Xgn q and Dw £ Xg^q-i- The argument 
for this is the same as for the n = 1 case [1]. 

The proof of uniqueness follows the same lines as in [1], by noting that the regularity 
requirements q > n/2 + 1 and qQ > n/2 + 1 + q are just sufficient in order to satisfy the 
hypothesis for q and qo in the fundamental energy estimate (cf. Lemma 3.3). □ 

3.3 Quasilinear symmetric hyperbolic Fuchsian systems 

We now have the tools needed to prove Theorem 2.4 for general quasilinear symmetric 
hyperbolic Fuchsian systems. The idea is to construct the following iteration scheme. 
We start with some seed function «r ] = no + u>[o] (here uq is the leading order term, and 
w\n\ is arbitrary) and we linearize the quasilinear system around u\py Then we apply our 
theory for linear equations above and hence, under the hypothesis of Theorem 2.4, this 
linear system has a solution um of the form um = uq + i%i . The next step is to linearize 
the quasilinear system around um. Again, we can apply the above theory to this system 
and hence obtain the solution up] = uo + topi- We thus obtain a sequence (t^ui) £ X$^ q 
with parameters 5, u, q as specified in the hypothesis of the theorem, which we seek to 
show converges to the remainder of the solution of the nonlinear equation. 

There are several steps to carry out in making this argument precise (and proving 
convergence). First, one constructs an operator G(uq) : B$^ qs — > Xg^g which maps 
iO[i] to W[i+i]- If the sequence (w[q) were to leave the closed ball Bs !fl , q , s , however, then 
only finitely many sequence elements would be defined. Using Eq. (3.6) and the first 
Lipschitz estimate in Condition (iii) of Theorem 2.4, it can be shown that this can be 
avoided. In fact, we can choose a sufficiently small 5' G (0, 6] so that G(ito) maps every 
w £ 13$/^^^ to an element in Bgi^ q s . Hence, if we restrict the elements Wu\ as above 
to the time interval (0,(5'], then the map G(zio) indeed generates a well-defined infinite 
sequence in B§i^ q%s . For every w\,W2 G Bg/^^g, we find a Lipschitz estimate of the 
form 

||G(« )[wi] - G(uo)[w2]\\s',n,q-l < C5' p \\wi - tf^lU'^g-l, 

for constants C, p > 0, thanks to the second estimate in Condition (iii) of Theorem 2.4. 
The reason for considering the (q — l)-norm as opposed to the (/-norm is that we are 
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required to compare a sequence of systems with different coefficients because the systems 
are obtained by linearizing around different solutions; recall that we are interested in 
quasilinear equations here. Just as in the proof of Proposition 3.5, this leads to a loss of 
one derivative. Therefore, since we apply Proposition 3.5 in each step of the iteration for 
solutions with q — 1 derivatives, we are led to a stricter inequality on q which is stated 
in the hypothesis of Theorem 2.4. We can assume without loss of generality that 5' is 
so small so that C5' p < 1. This 5' is referred to as 5 in the statement of Theorem 2.4. 
It follows that is a Cauchy sequence in Xs'^^-i with limit w. However, we know 

that the sequence (w^) is bounded in X^^^ q , and hence we can use duality arguments to 
show that w £ Xg/^^. It is then left to show that w is in fact differentiable in time and 
hence a solution of the equation. For this, we note that any solution to the nonlinear 
singular initial value problem is a fixed point of the operator G(uq). The limit w is the 
unique fixed point of G(uq) in Xs^^q, and hence it follows that w is the unique solution 
of the singular initial value problem. 

4 Concluding remarks 

In our previous work [1], we have used the results for one space dimension to construct 
a family of smooth (non-analytic) spacetime solutions to the 3 + 1 Einstein vacuum 
field equations with T 2 -symmetry, with these spacetimes all characterized by asymptot- 
ically velocity term dominated (AVTD) behavior in the neighborhood of the spacetime 
singularity. We expect that a similar conclusion should hold for [/(l)-symmetric space- 
time solutions and, in order to establish such a result, our generalization to Fuchsian 
systems with more than one spatial dimension should be crucial. As far as the ana- 
lytic class is concerned, earlier work by one of the authors with Moncrief and Choquet- 
Bruhat [9,10,13] uses the theorem in [17] to show that there are analytic [/(l)-symmetric 
solutions of Einstein's equations with AVTD behavior. 

We recall that a solution exhibits AVTD behavior if near the singularity the solu- 
tion evolves according to a simplified system, specifically an ordinary differential equa- 
tion (ODE), which can be derived from the full Einstein system by dropping spatial 
derivatives and keeping time derivative terms, only. The verification of such behavior 
has proven useful in establishing the validity of the strong cosmic censorship conjec- 
ture [11, 19], at least within restricted classes of solutions. We believe that this work 
could play a role in extending the class of solutions for which this conjecture is known 
to hold. 
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